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R. Isaacs, J. Ferrand, R. J. Duffin, and several others have developed a function theory for "discrete 
analytic functions" defined on "discrete regions" in the "discrete complex plane. " In this paper we 
bring to li^ht some combinatorial-topological properties of "simple discrete regions," and we study 
some basic properties of discrete analytic functions that are defined on simple discrete regions. These 
conihinatorial-topological properties and basic properties are then used to establish an existence and 
uniqueness theorem for discrete complex functions with prescribed "boundary values" and "residues" 
on an arbitrary simple discrete region. 
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1 . Introduction A y 

This paper is concerned with complex functions 
denned at the nodes of a square mesh Mh that lies 
in the complex plane and is depicted in figure 1.1. The 
theory of these functions was initiated by R. Isaacs 
[1] 1 and J. Ferrand [2 J and has been extensively de- 
veloped by R. J. Duffin [3] and several of his students 
[4,5,6,7]. The functions investigated are said to be 
"discrete analytic" and are discrete analogs of analytic 
functions of a continuous complex variable. A com- 
plex function / is discrete analytic at a square S be- 
longing to Mh if the difference quotient of / across 
one diagonal of S is equal to the difference quotient 
of /across the other diagonal of S. 

The results obtained in [1,2, and 3] include discrete 
analogs of the Cauchy Riemann equations, Laplace's 
equation, the maximum principle, conjugate harmonic 
functions, differentiation, integration, the residue 
theorem, Cauchy's theorem, Morera's theorem, 
Cauchy's integral formula, and polynomials. Some 
of the concepts studied have no direct analogy in 
the classical continuous theory. These include the 
notions of duality and bipolynomials. 

In [3 1 it is shown that for each real function u that 
is "discrete harmonic" on a "simple discrete region 
52" there is a real "conjugate function v that is 
discrete harmonic on 3& such that u+iv is discrete 
analytic on Si. This theorem becomes a corollary to 
the main theorem of this paper, which is stated and 
proved in section 3 as Theorem 3.5. The main theorem 
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Figure 1.1. The discrete complex plane. 

discusses both the existence and the uniqueness of 
a function u-\-iv that has prescribed "residue" and 
"boundary values" on a simple discrete region St. 

The proof of the main theorem depends on some 
properties of simple discrete regions and on some 
basic properties of discrete functions. The required 
properties of simple discrete regions are stated and 
proved in section 2, and the required basic properties 
of discrete functions are stated and proved in the first 
part of section 3. With the possible exception of Euler's 
formula [8] this paper is completely self-contained. 

2. The Discrete Complex Plane and 

Simple Discrete Regions 

The mesh Mh, as shown in figure 1.1, consists of a 
union of nonoverlapping squares, each of whose sides 
has length h. The lattice »§?/, (or discrete complex 
plane) is composed of the points that fall on a vertex 
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of a member of M n - Each member z of J^ i, can be rep- 
resented as z = xh-\~iyh where % and y are appropriate 
integers. The even lattice §T/, consists of those lattice 
points for which x + y is even, and the odd lattice €h 
consists of those lattice points for which x~\-y is odd. 

A chain is an oriented set of lattice points Zo, 
Zi, . . ., z« with the property that the straight line 
joining z\- \ to 2/ lies along an edge of a square belong- 
ing to JCu. A chain 2o, 2i, . . ., z„ for which zo — z n is 
called a simple closed chain when its elements, ex- 
clusive of 20 and 2> t , are distinct. An even chain (odd 
chain) is an oriented set of even (odd) lattice points 
2o, 21, . . ., z n with the property that the straight line 
joining z\-\ to z\ lies along a diagonal of a square 
belonging to M n - 

A set of lattice points is said to be connected if each 
pair of points belonging to the set can be connected 
by a chain so that each point of the chain lies in the 
set. A set of lattice points is said to be evenly con- 
nected (oddly connected), if each pair of even (odd) 
lattice points belonging to the set can be connected 
by an even (odd) chain so that each point of the even 
(odd) chain lies in the set. 

A discrete region in the discrete complex plane oSf/j 
consists of the nodes of a union of squares from Jth. 
The discrete region is said to be constructed from 
these squares, and the union of the closed planar sets 
bounded by these squares is called the associated 
region, The boundary of a discrete region con- 
sists of those lattice points that belong to the discrete 
region and lie on the boundary of its associated region. 
The union of all lattice points that belong to the dis- 
crete region but do not belong to its boundary is called 
the interior of the discrete region. Lemma 2.1 
is an obvious consequence of these definitions. 

Lemma 2.1: A lattice point z is an interior point of 
a discrete region & if, and only if, each of the four 
squares having z as a vertex belongs to the family of 
squares used in the construction of 01. 

A discrete region is said to be finite if it consists of 
a finite number of lattice points. A simple discrete 
region is a finite discrete region whose boundary can 
be represented by a simple closed chain. 

We shall now deduce some useful properties of 
simple discrete regions. It does not seem worthwhile 
to detail the combinatorial-topological machinery 
required for full formal rigor in the demonstration of 
these intuitively evident properties: instead, we aim 
at leaving the reader with a clear picture of just 
what 44 sub-lemmas" would have to be verified when 
using such machinery in a fully detailed treatment. 

Most of the proofs rely on induction. The following 
lemma provides the key to these proofs. 

Lemma 2.2: If & s is a simple discrete region con- 
structed from s squares and s ^2, then & s is the union 
of a simple discrete region & s -i constructed from 
s — 1 squares and a simple discrete region 3%\ con- 
structed from one square. 

PROOF: Let C: z , Zi, . . ., z n be the simple closed 
chain that represents the boundary of «^?. v , and suppose 
that C is oriented so that 3% % lies to the left as C is 
traversed from z<> to z n =Zo. Since C is finite and 
closed, there exists a lattice point z t eC: Zo, . . ., Zu 



. . . , z n at which C makes a left turn. The lattice can 
clearly be rotated so that Zt-i, Zt, and Z/+i have the 
positions shown in figure 2.1. According to our choice 
of C the square Si is contained in the associated 
region R s of Sft s . Since s 5= 2, at least one of the squares 
S», S3, and Si must be contained in R s . We must 
distinguish seven cases. 
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FIGURE 2.1. A left turn in a simple closed chain. 

Case I: S 2 , S3, and S4 are contained in R s . Then 
consider the closed chain C: 20, • • ">Zt^u zt, Zt+i 9 
. . ., z n . Since ~z t is an interior point of ^? s , C is simple 
and hence encloses a simple discrete region 9i s -\ 
C 0% 8 , In this case R\ = S\. 

Case II: S 2 and S3 are contained in R s but S4 is not. 
This means that Ti precedes zt-\ in C (Le., C: 20, . . ., 
Tu z t -i, z t , z t +i, . . ., z w ). Now let C: z , . . .,27, 
z*+i v . .,_z n be obtained from C by deleting Zt-\ and 
Zt. Then C is clearly simple and encloses a simple 
discrete region & s -\ C 3$ s . In this case R\=S\. 

Case III: S 2 and S 4 are contained in R s but S3 is not. 
This means that the sequence 6, Zt, d appears in C 
and hence C makes a right turn at zi. Since there must 
be four more left turns than right turns in C, we could 
have initially chosen zt so as to eliminate this case. 

Case IV: S3 and S 4 are contained in R s but S 2 is 
not. This case is handled in a manner similar to 
Case II. 

Case V: S2 is contained in R s , but S3 and S4 are not. 
This means that the sequence 6, 27, Zt-u %u Zt+i be- 
longs to C (i.e., C: 20 . . ., b, z t , Zt-i, z t , zt+\, . . ., z n ). 
Now let (7: 2 , . . . , />, z~t, Zt+i , . . . , z n be obtained 
from C by deleting zt-i and z t . Then C is clearly 
simple and encloses a simple discrete region 5? 5 _ 
C 8% 8 . In this case R\ = Si. 

CASE VI: S3 is contained in R s but S 2 and S4 are not. 
This case cannot occur, since Ft would have to appear 
twice in C, which would mean that C is not simple. 

Case VII: S 4 is contained in R s , but S 2 and S3 are 
not. This case is handled in a manner similar to Case V. 

All possibilities have been exhausted, so the proof 
of Lemma 2.2 is complete. 

We can now derive some useful properties of simple 
discrete regions. 

Lemma 2.3: If & s is a simple discrete region con- 
structed from s squares and s ^ 1, then & s is con- 
nected and both evenly and oddly connected. 
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PROOF: The theorem is obviously true for 5=1. 
If 5 > 1, then Lemma 2.2. shows that & s = & s -i U 9i u 
Using induction, we assume that both 9$ 8 -\ and £%\ 
have the desired properties. Since 3t 8 is simple, 
3fts-\ and £fc\ must have at least one even lattice point 
in common and at least one odd lattice point in com- 
mon. The rest of the proof is trivial. 

Lemma 2.4: // «^? s is a simple discrete region con- 
structed from s squares and Si is one of the squares 
used in constructing & s , then Si is the first member 
of a sequence {Sj}*j = 1 of squares belonging to M^ 
such that 

s 

(1) «^ s is constructed from U Sj, and 

j=i 

(2) For each j such that 1 < j ^ s, there is an integer 
k such that 1 ^ k < j and such that Sj has a side 
in common with S k . 

PROOF: The proof, being trivial for 5=1 and 5 = 2, 
is accomplished through induction on 5. If 5^=2, 
then Lemma 2.2 shows that ^, = ^,_, U 9t x . Two 
cases must be treated. 

Case I: Si is used in the construction of ffl s -\. Then 
there exists a sequence {Sj}j'= / that satisfies condition 
(2) and is used for the construction of 9t % ~\. Letting 
S s be the square used in constructing &\ and re- 
membering that the boundary of 9i s is simple, we see 
that Si must have a side in common with a member of 
{SjYjz}, which means that {Sj}^^ , satisfies both 
condition (1) and condition (2). 

Case II: Si is used in the construction of 3i\. Since 
the boundary of 3%^ is simple, we can choose a square 
S-> that is used in the construction of 3$ s -\ and has a 
side in common with Si. Then there is a sequence 
{Sj + \}j=l that satisfies condition (2) and is used for 
the construction of {% 8 - \. It is now easily seen that the 
sequence {Sj}j_ ! satisfies conditions (1) and (2). 
This completes the proof of Lemma 2.4. 

Lemma 2.5: Let p be the number of lattice points 
in a simple discrete region & s that is constructed 
from s squares and has b boundary points. Then 
2p - b - 2 = 2s. 

PROOF: The proof of this lemma can be accom- 
plished through induction on 5. A shorter proof can 
be based on Euler's formula from the theory of planar 
graphs [8j. If a planar graph is connected, the Euler 
formula states that 



e+f=2, 



(2.1) 



where v is the number of vertices in the graph, e is 
the number of edges in the graph, and /is the number 
of regions into which the plane is divided by the 
graph. In order to apply this formula we must con- 
sider the connected planar graph G that results from 
taking the union of the squares used in the construc- 
tion of 3$ s . The number of vertices in G is simply the 
number of lattice points in 91^ that is, 



Since 5 squares are used in constructing G, the number 
of regions into which the plane is divided is 5+1; 
thus 



/=* + !. 



(2.3) 



There are b edges connecting the boundary points of 
5? s , because the boundary of 3$ s is a simple closed 
chain. Each of these edges borders only one of the 
squares used in the construction of @% 8 . All other edges 
belonging to G border two of the squares used in the 
construction of 3£ 8 . Since each square is bordered by 
4 edges, it follows quite easily that 



2e = 45 + b. 



(2.4) 



Multiplying eq (2.1) by 2 and substituting eqs (2.2), 
(2.3), and (2.4) into the result proves Lemma 2.5. 

Figure 2.2 shows the union of ten squares to form 
a simple discrete region ^?io. Odd lattice points are 
denoted by crosses, and even lattice points are rep- 
resented by dots. The square labeled with Sj occupies 
the /th position in the ordered sequence {Sj}j ( l iof 
squares used to construct & ]{) . The points z\, z 2 , z.?, 
and za form a chain. An even chain zf, z$, zf, z e A is 
shown along with an odd chain z\\ z§, z%. The boundary 
of ^io, denoted by d$?io, is depicted by the simple 
closed chain z u z 2 ,. . ., Zie- The node z\ is an example 
of an interior point of &\ {) . 
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FIGURE 2.2. Chains and a simple discrete region. 

3. The Main Theorem 

We begin by defining a finite difference operator 
L whose definition is taken from [3]. 

Def. 3.1: Let f be a complex lattice function and 
suppose that S is a square belonging to Ji h . Then the 
residue L(f, S) of f at the square S is defined by 

L(f, S) = f(z ) + ifizi) + i 2 f(z 2 ) + i 3 f(z,), 



107 



where zo, Zi, Z2, a/id z 3 /w* f/ie orientation depicted 
in figure 3.1 and i = V— 1 . 
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Figure 3.2. An interior point. 



Figure 3.1. Orientation of the square S. 

The quantity L(f,S) is termed the residue of / at S, 
because in "discrete contour integration" it plays a 
role [3] analogous to the role played by the residue of 
a complex function of a continuous complex variable 
in the classical theory of contour integration. 

Def. 3.2: A lattice function f is discrete analytic 
at a square S belonging to Jt h if its residue 
L(f, S) at S is zero. A lattice function f is discrete 
analytic in a discrete region if it is discrete 
analytic at each square used in the construction of ' 0t. 
The reader should have no trouble showing that this 
definition of discrete analyticity is equivalent to the 
one given in the first paragraph of section 1. 

To find analogues of the Cauchy Riemann equations, 
one need only set the real and imaginary parts of 
L(f,S) equal to zero. This gives the following theorem, 
which is taken from [3]. 

THEOREM 3.1: Let f=u + iv where u and v are real 
lattice functions. Then f is discrete analytic at SeJ? h 
if and only if 



(i) 
and 



(ii) 



u(z 2 ) — u(z ) = v(z 3 ) — v(zi 



v(z 2 ) — v(z ) = — [u(z 3 ) - u{zi) 



where z , Zi, Z2, and z 3 have the orientation shown by 
square S of figure 3.1. 

Following the example of the theory of functions of 
a continuous complex variable, we give a definition 
of discrete harmonicity. This definition, which is 
taken from [3], is justified by Lemma 2.1. 

Def. 3.3: Suppose f is a complex lattice function 
defined on a discrete region 01 of the discrete complex 
plane J?f h . Let a function H(f, •), which maps the in- 
terior points of into the field of complex numbers, 
be defined by 

H(f, Zi ) ^ f(zi) + f(z 2 ) + f(z 3 ) + f(z 4 ) - 4f( Zi ) 

where z i? zi, z 2 , z 3 , and z 4 have the positions shown in 
figure 3.2. Then f is said to be discrete harmonic 
at an interior point z x of if H(f, Zj) = 0. 

A simple consequence of this definition is the "max- 
imum principle" given by Theorem 3.2. 



THEOREM 3.2: Ifiis a complex lattice function that 
is discrete harmonic at each interior point of a simple 
discrete region 0, then |f | assumes its absolute max- 
imum on the boundary of 0. 

PROOF: Since 01 is finite, |/| assumes an absolute 
maximum Mat a point z\ of 0. If Z\ is on the boundary 
of ^?, the theorem is proved. 

Suppose then that z\ is in the interior of ffl. If Z\ is 
even (odd) then Lemma 2.3 enables us to construct 
an even (odd) chain C: zi, 2:2, . . ., z n from z\ to an 
even (odd) boundary point z n of 01 with the property 
that Zi, f = l, 2, . . ., n — 1 is an interior point of 01. 
Let j be the smallest positive integer less than or 
equal to n for which \f(z;)\ < M. Sincey" must be greater 
than one, we know that zy-\ exists and is an interior 
point of 01 with the property that |/Uj_i)| =M. The 
lattice J? h can clearly be rotated so that Zj-\ and Zj 
have the orientation shown in figure 3.3. 



Z j Z 3 



Figure 3.3. An interior point on an even (odd) chain. 



Definition 3.3 shows that 



/te-i 



_ /(zj)+/(z.)+/fa)+/(2») 
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Hence 



\f(zj-i)\ 



[f(zj)\ + [flz,)| + [f(z-M 
4 



[f(zj)\<[f(zj- l )\=M, 



\f(zs)] 
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we see that 



f( Zl ) +/(z 3 ) +/(z 5 ) +/Ut) -4/(z,) =0. q.e.d. 



M < = M, 



which is a contradiction. Hence [f(zj)\=M forj=l, 
. . ., n and thus \f(z n )\ — M. This completes the proof 
of Theorem 3.2. 

We can use the preceding theorem to discuss the 
"Dirichlet problem" for the lattice J2V 

THEOREM 3.3: Suppose 0Z is a simple discrete region 
containing i interior points along with b boundary 
points. //{Bi, B2, . . ., B b } is a set ofb complex num- 
bers, then there is a unique lattice function f that is 
discrete harmonic in St and assumes the values Bi, 
B2, . . ., B b on the boundary of St. 

PROOF: Definition 3.3 provides us with a system of 
i linear equations in i unknowns, which must be 
satisfied by a discrete harmonic function defined on St. 
The corresponding homogeneous system is obtained 
by setting each boundary value of / equal to zero. 
Theorem 3.2 shows that the homogeneous system has 
only the trivial solution. Hence, the inhomogeneous 
system has a unique solution, which proves the 
theorem. 

The next theorem, which is taken from [3], relates 
the concept of discrete harmonicity to the concept 
of discrete analyticity. 

Theorem 3.4: If f is a lattice function that is dis- 
crete analytic at each square making up a discrete 
region St, then f along with its real and imaginary 
parts is discrete harmonic at each interior point of SI. 

PROOF: Suppose zj is an interior point of S. Then, 
according to Lemma 2.1, the squares Si, S2, S3, and 
Sa shown in figure 3.4 belong to the family of squares 
used in the construction of 3i. Hence 

f( Zl ) + if(z 2 ) + iY(zj) + fifos) = 
f{z 2 ) + if (z») + ij{z 4 ) + ij(zj) = 
f(zj) + if(zd + ij{z- a ) + i*f(ze) = 
f(zs) + if(zj) + i*f(z«) + Pf(z 7 ) = 0. 

Multiplying the first equation by + 1, the second equa- 
tion by — i, the third equation by —1, and the fourth 
equation by + i, and then adding the results shows that 
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The preceding theory can be used to discuss the 
existence and uniqueness of a lattice function / that 
is defined on a simple discrete region St and has pre- 
scribed residue at each square making up St. The 
following existence and uniqueness theorem is the 
main result of this paper. 

THEOREM 3.5: Let St be a simple discrete region 
constructed from s squares, and suppose that St has 
p lattice points and b boundary points. If 

(i) z° and z e are odd and even lattice points respec- 
tively belonging to the discrete region St, 
(ii) {v°, v e } is a set of two real numbers, 

(iii) {u J }j 3 =] is a set ofb real numbers, and 

(iv) {a J + ib J }j = 1 is a set of s complex numbers, 
then there is a unique complex function f=u + iv that, 
is defined on St and has the following properties: 

(I) v assumes the values v° and v e at z° and z e , 

respectively. 
(II) u assumes the value u j at the jth boundary point 

of St. 
(Ill) L(f, S J ) = a J +ib J /orj = l, . . ., s where {S J }jL, 

is the set of squares used in the construction 

of St. 

PROOF: Equating the real and imaginary parts of 
the equation L(f, &) = at ' + ib* shows thai f= u + iv 
has residue a j -f ib> at & if, and only if, 



and 



u(z 2 ) — u(zo) =v(z :i ) -v(zi) —a) 



v{z 2 ) —v(zo) =— [u(z:i) —u(zi) +#], 



where zo, Zi, z 2 , and 23 have the orientation shown in 
figure 3.1. Since St is constructed from s squares, we 
obtain a system of 25 linear equations that must be 
satisfied by the real and imaginary parts of the func- 
tion / whose prescribed residues are {a) + ibf)^ v 
There are 2p — b~2 unknowns because u is given at 
b boundary points and v is specified at z° and z € . 
Lemma 2.5 shows that the number of unknowns is 
equal to the number of equations in the system. The 
corresponding homogeneous system is obtained by 



setting u l , u z , . 



u u , v v , v*, a 1 



., a s , 6 1 , . . ., 



FIGURE 3.4. An interior point of a discrete region. 



b s ~ l , and b s equal to zero. 

Suppose u and v represent a solution of the homo- 
geneous system. Then/=a-h/t> is discrete analytic 
on St, and Theorem 3.4 shows that u is discrete har- 
monic at interior points of St. We conclude from The- 
orem 3.3 that u is identically zero. Now suppose that 
z° is a vertex of a square Si and that {Sj}j =1 is the cor- 
responding sequence of squares given by Lemma 2.4. 
Since u is identically zero and v(z°) =0, the "Cauchy 
Riemann equations" of Theorem 3.1 show that v is 
zero at the other odd lattice point of Si. Using the 
second conclusion of Lemma 2.4, we observe that 
.v is, a fortiori, zero at one of the odd vertices of S2 and 
hence zero at the other odd vertex of S 2 .- Repeating 
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this argument an additional (5 — 2) times, we conclude 
that v is zero at all odd lattice points of &. The same 
argument with z° replaced by z e and the word odd re- 
placed by the word even shows that v is zero at all 
even lattice points of 8&. Hence, both u and v are 
identically zero on 01 and thus the homogeneous sys- 
tem has only the trivial solution. This means the inho- 
mogeneous system has a unique solution, which 
proves the theorem. 

COROLLARY 3.1: If u is a real discrete harmonic 
function that is defined on a simple discrete region ^?, 
then u is the real part of a lattice function u + iv that 
is discrete analytic on Si. 

PROOF: Theorem 3.5 shows the existence of a dis- 
crete analytic function u + iv with u assuming the 
boundary values taken on by u. Then Theorem 3.3 

shows that u = u. 

The real function v of Corollary 3.1 is called a dis- 
crete harmonic conjugate of u. In the theory of 
functions of a continuous complex variable a particu- 
lar harmonic conjugate of a real harmonic function is 
uniquely determined by its value at a given point of 
its domain. An inspection of Theorem 3.5 shows that 
a particular discrete harmonic conjugate v of a dis- 
crete harmonic function u is uniquely determined by 
its values at two different points z° and z e of its domain. 



The difference between the two theories is due to the 
fact that the discrete harmonic operator H does not 
relate a discrete function restricted to the odd lattice 
to the discrete function restricted to the even lattice. 
In solving the Dirichlet problem for J?/,, we must, in 
a sense, solve two Dirichlet problems, one for the 
even lattice and one for the odd lattice. 
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